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ABSTRACT 

Ultraviolet observations of classical T Tauri Stars (cTTSs) have shown that there 
is a hot (T ~ 80, 000 K) and dense (n e ~ 10 10 cm -3 ) component associated with 
the large scale jet. This hot component is formed very close to the base of the jet 
providing fundamental information on the jet formation mechanism. In this series, 
we have investigated whether this component can be formed in disc winds, either 
cool or warm. To conclude the series, jet launching from the interface between the 
magnetic rotor (the star) and the disc is studied. Synthetic profiles are calculated 
from numerical simulations of outflow launching by star-disc interaction. Profiles are 
calculated for several possible configurations of the stellar field: dipolar (with surface 
strengths, £>» of 1, 2 and 5 kG) or dynamo fed. Also two types of discs, passive or 
subjected to an af2-dynamo, are considered. These profiles have been used to define the 
locus of the various models in the observational diagram: dispersion versus centroid, 
for the profiles of the Si III] line. Bulk motions produce an increasing broadening of 
the profile as the lever arm launching the jet becomes more efficient; predicted profiles 
are however, sensitive to the disc inclination. Models are compared with observations 
of the Si III] lines obtained with the Hubble Space Telescope. 

In addition, it is shown that the non-stationary nature of star-disc winds produce a 
flickering of the profile during quiescence with variations in the line flux of about 10%. 
At outburst, accretion signatures appear in the profiles together with an enhancement 
of the wind, producing the correlation between accretion and outflow as reported from 
RU Lup, A A Tau and RW Aur observations. 

Key words: stars: pre-main-sequence - stars: winds, outflows - stars: circumstellar 
matter - stars: formation 



1 INTRODUCTION 

Star formation is accompanied by heavy mass outflow. Out- 
flows are powered by the gravitational field (through the 
centrifugal gear) and by the release of the magnetic energy 
built up in the sheared and turbulent protostellar discs. 
Since the very early studies of the atmospheres of the T 
Tauri Stars (TTSs), it is known that during the classical 
T Tauri phase (~ 1 Myr old stars), TTSs drive powerful 
and massive cool (T~ 10, 000 K) winds whose signature was 
first detected in the profiles of neutral and singly ionized 
species being the most conspicuous the optical Ha line and 
the Mg II resonance multiplet (UV1) in the UV (Giampapa 
et al 1981, Penston & Lago 1983, Finkenzeller & Mundt 
1984, Calvet et al 1985). The realisation that TTSs drive 
cool (T~ 10, 000 K) jets and bipolar outflows came in the 
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early 80's (i.e. Schwartz 1983, Mundt et al 1987). Ever since, 
it is unclear which fraction of cool wind is associated with 
the so-called cool stellar wind (and which fraction is associ- 
ated with the bipolar outflow). For instance, the absorption 
by this cool wind of the FU emission from the disc (produced 
within ~ 2 AU from the star, Herczeg et al 2004), clearly 
indicates that a significant fraction of the absorbing gas is 
associated with the large scale jet. 

Though there is a general consensus on bipolar out- 
flows being driven by the disc-star interaction and the disc 
wind (see i.e. Ferreira et al 2006), the nature of TTSs stel- 
lar winds remains open. Hot winds are expected to be pro- 
duced by the star since the strong magnetic activity and the 
powerful stellar corona (Preibish et al 2004) suggest that a 
pre-main sequence analogue to the solar coronal winds must 
exist. Unfortunately, this hot component has not been con- 
clusively detected. This is caused by the complexity of the 
atmospheres and the disc-star interaction region in TTSs 
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that makes difficult to disentangle the contribution from the 
photoionized accretion flow and the inner disc from a theo- 
retically hypothesised hot stellar wind. In 2005, Dupree et 
al, claimed to have detected, for the first time, this hot wind. 
These authors interpreted the absence of blueshifted emis- 
sion in the O VI resonance lines in TW Hya and T Tau, as an 
indication of the existence of a hot wind with temperatures 
as high as 300,000 K and mass-loss rates of 2 x 10 _11 M Q yr~ 1 
that absorbs the blue wing of the profile. As pointed out by 
Johns-Krull & Herczeg (2007), the lack of a blueshifted com- 
ponent could just be an indication of the line being formed 
in the accretion flow. Matter infall onto the stellar surface 
is channelled by the magnetic field driving to the formation 
of shocks at the impact points where the kinetic energy of 
the infalling gas is finally released. The temperature reached 
at the shocks can be as high as some 10 6 K and the pho- 
toionizing X-ray radiation preionizes the infalling gas col- 
umn that radiates over a wide range of temperatures and 
tracers (Lamzin 1998). 

An independent line of research, the one addressed in 
this series, was to investigate whether the base of the jet 
could be hot enough to radiate at temperatures as high as 
50,000 K-80,000 K, well below the temperature range of the 
hot stellar wind but above the fiducial temperatures of opti- 
cal jets: 10,000 K-20,000 K. The motivation for this research 
came from the discovery of C III] and Si III] semiforbidden 
emission in the UV spectrum of RY Tau and RU Lup at the 
same velocity of the optical jet; this emission traces gas at 
logT e ~ 4.6-4.8 (Gomez de Castro & Verdugo 2001). Cen- 
trifugally driven flows from magnetised accretion discs are 
submitted to pinching stresses j x B because the toroidal 
magnetic field (B) collimates the current of gas along the 
disc axis (j). Recollimation can drive to the formation of 
focal surfaces or shocks on the jet axis that are able to heat 
the gas to the some 10,000 K traced by the optical forbidden 
lines of S II] (Gomez de Castro & Pudritz 1993). However, 
temperatures as high as 80,000 K cannot be produced in 
cool disc winds. In the first article of this series, we ex- 
amined whether the photoionisation of cool disc winds by 
the stellar corona could cause the observed emission (Fcrro- 
Fontan & Gomez de Castro 2003, hereafter Paper I). We 
found that the propagation of the stellar radiation generates 
a cocoon of photoionized gas around the star. The extent of 
the photoionized region is small (tenths of au) in dense out- 
flows and close to the disc plane; however, it may cover the 
whole wind extent in diffuse winds, e.g. disc winds generated 
by small accretion rates (^ 10 -9 MQyr" 1 ). Photoionisation 
also modifies the electron density in the plasma. The in- 
terplay between ambipolar diffusion and the radiation field 
controls the electron temperature of the wind that is kept 
around 10,000 K, well below the temperatures traced by the 
UV semiforbidden lines. 

Disc winds are a fundamental mechanism for angular 
momentum transport in protostellar discs, contributing to 
the regulation of the accretion rate onto the star. However, 
most of the transport occurs in the disc itself through the 
magnetorotational instability (Balbus & Hawley 1991): weak 
fields provide a tension force that allows two orbiting fluid el- 
ements to exchange angular momentum on larger scales than 
the hydrodynamical viscosity scales. The magneto rotational 
instability acts like a dynamo, amplifying the field which is 
lost due to buoyancy, leading to a magnetised corona. Flares 



associated with reconnection events would naturally be pro- 
duced leading to the formation of a hot disc corona. Warmer 
disc winds, centrifugally launched from the disc corona, were 
shown to be able to reproduce the observed line fluxes and 
line ratios provided that the winds are clumpy (with filling 
factors about 1%) in Gomez de Castro & Ferro-Fontan 2005 
(hereafter Paper II), the second article of this series. Warm 
disc winds have also been shown to be able to reproduce 
larger scale jet observations (see i.e. Vlahakis & Tsinganos 
1999, Ferreira 2004). Warm disc winds provide an elegant 
and simple solution to the high jet temperatures observed 
in the UV however, they are not adequate to reproduce in 
full detail the observed physics of the line formation region. 
Long tails of bluewards shifted emission are detected in the 
Si III] and C III] plasma tracers, suggesting that the outflow 
launching mechanism is more efficient and, at the same time, 
the outflow is less collimated than predicted by the simple 
self-similar warm disc wind theory. As shown in Paper II, 
the self-similar solutions that are able to reproduce the ob- 
served jet properties (velocity and mass flow) require that 
the ratio between the sound speed and the escape velocity is 
0.43 at the Alfven radius. As a result, the radial expansion of 
the outflow - the main source of line broadening - is shifted 
to distances of about 10 AU from the star, where plasma 
cannot contribute to the Si III] and C III] emission because 
the densities and temperatures are very low and radiative 
cooling is dominated by singly ionized species (see Fig. 3 
in Gomez de Castro & Verdugo 2007, hereafter GdCV07). 
Thus, self-similar models produce winds too collimated at 
the base to reproduce the observed broadening of the Si III] 
and C III] profiles. 

The high densities revealed by the Si III], C III] and 
C IV] line ratios indicate that this line radiation is produced 
very close to the star. Henceforth the radiating plasma must 
be strongly affected by the disc-star interaction. In fact, 
there is evidence that Si III] and C III] radiation can be 
produced in structures very close to the star such as the 
ionized plasma torus around RW Aur (Gomez de Castro & 
Verdugo 2003) or the accreting shell in RY Tau (GdCV07). 
The physics of the interaction between the stellar magnetic 
field and the accretion disc is very rich, and it has been 
shown extensively that jet launching can be produced from 
this interaction region with physical temperatures and densi- 
ties similar to those of stellar chromospheres (Goodson et al. 
1997, von Rekowski & Brandenburg 2004, 2006). Jet launch- 
ing from the interface between the magnetic rotor (the star) 
and the disc is much more efficient than a pure disc wind, 
firstly because the centrifugal gear is higher closer to the star 
and secondly because of the heavy mass load onto the field 
lines at this point. Thus, to conclude this series, we have 
computed the radiative output from numerical simulations 
of outflow launching from star- disc interaction. It will be 
shown that the radial flow expansion, i.e. the outflow from 
the disc-star interaction region, is able to reproduce the ob- 
served broadenings for magnetospheric fields of about 1-2 
kG. In Section 2, the characteristics of the simulations from 
von Rekowski & Brandenburg (2004, 2006, hereafter vRB04 
and vRB06, respectively), that are used in this work, are 
summarised. In Section 3, the procedure to derive line fluxes 
and profiles from the simulations is described. The analysis 
of the output profiles is developed in Section 4. A discussion 



Dense, hot winds from TTSs 3 



on the implications and relevance of the results is deferred 
to Section 5, as well as the main conclusions. 



2 THE T TAURI STAR-DISC SYSTEM 

2.1 The magnetospheric star-disc model 

The interaction between the star and the Inner disc Bound- 
ary LayeiQ generates a sheared region fed by turbulent, mag- 
netised material from the accretion disc. Shear amplifies the 
stellar magnetic field, producing a strong toroidal magnetic 
field component. This toroidal field and the associated mag- 
netic pressure push the stellar poloidal field away from the 
stellar/disc rotation axis, inflating and opening the poloidal 
field lines in a butterfly like pattern thus producing a current 
layer between the stellar and the disc dominated regions as 
displayed in Fig. I. The magnetic link between the star and 
the disc is broken and re-established periodically by mag- 
netic reconnection. The opening angle of the current layer, 
as well as its extent, depends on the stellar and disc fields, 
the accretion rate and the ratio between the inner disc radius 
and the stellar rotation frequency. 

In the magnetospheric models of vRB04, a slow, hot 
and dense outflow driven mostly by poloidal magnetic pres- 
sure is emanating from stellar regions close to the rotation 
axis while fast, cooler and less dense magneto-centrifugally 
accelerated outflows are emanating from lower stellar lati- 
tudes and from the inner disc. The main difference with the 
magnetospheric models of Goodson et al. (1997) is that in 
the Goodson models the axial (stellar) jet is fast and well- 
collimated - and driven by magneto-centrifugal processes 
(see also Hirose et al. 1997) - and the inner disc wind is 
divergent, whereas in vRB04 the stellar wind is slow but 
on the other hand some collimation is clearly seen in the 
inner (magneto-centrifugally accelerated) disc wind. How- 
ever, the stellar wind seen in Romanova et al. (2002) is also 
slow, travelling into a rarefied corona. A very fast (warm 
and dense) stellar wind is seen in stellar dynamo models of 
vRB06; this wind is mostly driven by toroidal and poloidal 
magnetic pressure, as well as gas pressure. 

2.2 Basic characteristics of the numerical models 

Though a full description of the numerical simulations can 
be found in vRB04 and vRB06, we should briefly summarise 
their main properties. In both papers, the evolution of the 
flow, magnetic field, density and temperature is found by 
solving the continuity, the Navier-Stokes and the mean- 
field induction MHD equations for an axisymmetric system 
in cylindrical polar coordinates, assuming a piecewise poly- 
tropic model. Dynamo action in the disc (if present) is pre- 
scribed by a standard a 2 Q dynamo (e.g. Krause & Raedler 
1980) , where a is the mean-field a-effect and Q is the angu- 
lar velocity of the orbiting gas. a-quenching is included so 
that the disc dynamo saturates at a level close to equiparti- 
tion between magnetic and thermal energies. The code uses 

1 The shearing layer between the region dominated by the stellar 
magnetic field and the gravitationally dominated Keplerian disc 
has a finite extent. 
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Figure 1. The interaction between the stellar magnetic field and 
the disc twists the stellar field lines due to the differential rota- 
tion. The toroidal magnetic field generated out of the poloidal flux 
and the associated magnetic pressure tends to push the poloidal 
field lines outwards, inflating them, and eventually braking the 
magnetic link between the star and the disc (boundary between 
regions I and II). Three basic regions can be defined: Region I 
dominated by the stellar wind, Region II dominated by the disc 
wind and Region III dominated by stellar magnetospheric phe- 
nomena. The dashed lines trace the boundaries between these 
three regions. The continuous lines indicate the topology of the 
magnetic field and the shadowed areas represent regions where 
magnetic reconnection events are likely to occur, producing high 
energy radiation and particles (from Gomez de Castro 2004). 



dimensionless variables that have been scaled using as ref- 
erence values a typical sound speed of the coronal gas (100 
km s _1 ) and a typical surface density at the surface of the 
disc (1 g cm -2 ). Furthermore, the mass of the star has been 
taken as 1 Mq (solar mass), the mean specific weight as 
H — 0.6 and the polytropic index as 7 = 5/3. The compu- 
tations have been carried out in a domain of extent 0.2 AU 
in the radial (ro) direction and ±0.1 AU above/below the 
disc mid-plane; mesh sizes are 5zu — 8z — 0.001AU. The 
inner edge of the disc is at 4 stellar radii (in the models 
with stellar dipolar magnetosphere; cf. Table 1) or at 2.4 
stellar radii (in the model with stellar dynamo), which in 
both cases corresponds to 12 solar radii. The disc extends 
to the outer boundary of the computational domain. 

The vRB04 models represent a step forward over the 
self-similar warm disc winds models of Paper II. Around the 
inner edge of the disc, the interaction with the stellar mag- 
netosphere and the stellar wind is taken into account. Also, 
the rigidity of the self-similarity constraint is lost. Further- 
more, the Lorentz force j X B is included in the Navier- 
Stokes equation, as well as the generation of magnetic fields 
by the standard o?Q, dynamo in the disc. The saturation 
level of the disc dynamo is governed by the equipartition 
between magnetic and thermal energies in the disc. In ac- 
cordance with observations of protostellar star-disc systems, 
the model implements a dense, cool disc embedded in a rar- 
efied, hot disc corona. This is achieved by prescribing an en- 
tropy contrast between the disc and corona so that specific 
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entropy is smaller within the disc and larger in the corona, 
and by choosing hydrostatic equilibrium as the initial state. 
Therefore, detailed modelling of the coronal heating physics 
is avoided by prescribing a "thermal" frontier between disc 
and corona that is based on the observational properties of 
accretion discs (piecewise adiabatic model in disc, corona 
and star). However, temperature is evolving with time, in 
particular in the disc corona as one can clearly see in the 
relevant figures. Assuming a perfect gas, the continuity equa- 
tion is re-formulated as an evolution equation for specific en- 
thalpy, which is directly related to temperature. The specific 
enthalpy is then combined with the gravitational potential 
to form the potential enthalpy. 

In a sense, the disc in the vRB04 simulations can be 
understood as a dynamical boundary that regulates the ex- 
tra electromotive force caused by the disc dynamo and taps 
the mass flow into the wind and onto the accreting star. 
vRB04 adds to the dynamics of the self-similar models in 
Papers I & II but shares a similar thermal treatment allow- 
ing a straightforward comparison. 

The full set of simulations in vRB04 and vRB06 ex- 
plores the role of disc-star interaction in mass ejection for 
various configurations of the stellar field (a dipolar mag- 
netosphere or a mean-field a Q dynamo generated field) 
and strengths (models with different stellar surface fields 
are analysed). The simulations also explore the role of mag- 
netically active discs (with a 2 Q dynamos) in the dynamical 
evolution of the system by comparing the active disc evo- 
lution with that of passive discs. In all the simulations, the 
three basic winds: hot stellar wind, warm disc wind and 
episodic ejections from reconnection in the star-disc current 
layer are readily identified. In this work the radiative output 
has been calculated only for a subset of these simulations. 

Firstly, we have selected a reference model to analyse 
whether and how the temporal evolution of the simulated 
star-disc system shows in the profiles of the spectral tracers 
of the outflows. This is model "Ml" in vRB04. The stel- 
lar magnetosphere is initially a dipole threading the disc 
but with time can freely evolve outside the anchoring re- 
gion away from the stellar surface; the dipole is aligned with 
the disc axis and has a strength at the stellar surface of 
B, ~ lkG. The disc is assumed to be undergoing a signif- 
icant dynamo action with ao = —0.1. This is the so-called 
"Reference" model in Table 1. The thermal properties of the 
disc wind at the launching points in the corona are similar 
to those derived in Paper II for warm disc winds, i.e.: 

T = ui^H^)©-' m 

where M* is the mass of the T Tauri star, u>o is the radial 
distance of the wind foot point (in the disc corona) and y 
is the ratio between the outflow density at vjo and at the 
Alfven radius (see Paper II for more details). The temper- 
ature drops as vj^ 1 at the base of the disc wind, similarly 
to what derived in vRB04 for the Reference model (see i.e. 
Fig. 2 in vRB04). Notice that in Paper II these conditions 
are set to reproduce the observed large scale properties of 
protostellar jets following Vlahakis & Tsinganos (1999). In 
vRB04, the thermal properties of the wind are derived from 
the specific enthalpy, h, that is included in the modified 
Navier-Stokes and continuity equations (von Rekowski et 
al. 2003). The system is piecewise adiabatic and the equa- 



Table 1. Summary of the simulations analysed 



Iden 


Stellar field 




disc field 




Reference 


Dipolar, B, <- 


- lkG 


a 2 f2 dynamo with ao = 


-0.1 


Mag-2kG 


Dipolar, B* /• 


- 2kG 


a 2 f2 dynamo with ao = 


-0.1 


Mag-5kG 


Dipolar, B* <- 


- 5kG 


a 2 f2 dynamo with ao = 


-0.1 


Passive 


Dipolar, B, r 


- lkG 


No dynamo 




Dynamo * 


a 2 f2 dynamo, ao = 1-5 


ct 2 £l dynamo with ao = 


-0.15 



tion of state a polytrope, thus the sound speed is given as 
c 2 = (7— l)h, with 7 = 5/3, and the temperature is derived 
directly from h = c p T. 

In addition, we have selected a set of models to study 
the variation of the quiescent profiles with the strength of the 
magnetospheric field. Models with B* ~ 2kG and B, ~ 5kG 
have been analysed. These models are quoted as "Mag-2kG" 
and "MagSkG", respectively. As in the reference model, the 
disc is assumed to be undergoing dynamo action with ao = 
-0.1. 

The next step has been to explore whether the presence 
of an active or passive disc can be detected in the profiles. 
For this purpose we have analysed a simulation similar to 
that of the "Reference" model but with a passive disc (no 
disc dynamo); this model is quoted as "Passive" in Table 1. 

Finally, we have analysed how the profiles change with 
the characteristics of the stellar field, i.e. we have inves- 
tigated whether introducing an active dynamo in the star 
modifies the line profiles. This model is quoted as "Dy- 
namo*"; the star is considered to be fully convective with 
the stellar dynamo rooted up to about 5 solar radii and 
ao = +1.5 for the stellar dynamo. 

A summary with the names and main magnetic prop- 
erties of the models is given in Table 1. 

The vRB04 models analyse the dynamical evolution 
of the star-disc interaction during total run times between 
~ 150 d and ~ 1900 days. However, the time scale anal- 
ysed in this work is significantly shorter: about 75 d in 
the "Reference" model where the effect of outburst com- 
pared to quiescence is studied in the profiles. This shorter 
time scale is, in general, about twice the typical matter 
crossing time from the outer to the inner disc border: 
(t*7 ou t — Win )/u zu «(0.2 - 0.06)AU/7kms- 1 ~ 34.5 d. In 
vRB04, matter is injected locally wherever and whenever 
it is needed in the disc region, in order to model mass ac- 
cretion from the outer parts of the disc to the inner parts. 
This might produce artificially high mass loss rates in the 
long term, especially around the inner disc edge where mass 
outflow is more efficient (see vRB04 for a detailed discus- 
sion on the implications of this approach). Selecting periods 
shorter than the crossing time minimises the effect of matter 
injection. 



3 COMPUTATION OF THE LINE PROFILES 

Flow kinematics is handy described by three components: 

• Motion along the rotation (or 2) axis (vertical flow com- 
ponent): V z 

• Rotation given by the toroidal flow component: Vt 
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Figure 2. Si III] emissivity maps of the reference model during 
quiescence (top) and in outburst (bottom). Notice that at out- 
burst the Si III] emitting region is large and reaches closer to the 
star. 



• Expansion away from the rotation axis (radial flow com- 
ponent): V<sj 

The numerical simulations provide as output the density pk, 
temperature Tk, the three velocity components (14,fc> K,k, 
V roj ft) and the three magnetic field components at any given 
cell, k, of the computational mesh. As the models are ax- 
isymmetric, the simulations provide the same temperature 
and density for a cylindrical ring around the disc axis. As 
the C III] 1908 and the Si III]i892 lines are semiforbidden, the 
output radiation can be directly computed from the emis- 
sivity and the volume of the ring. The spectral distribution 
of the radiation is then computed for a given fc-ring using 
the projection of the ring velocity into the line of sight. At 



oint in the ring, the velocity projected into the line of 
is given by 

= —Vm,k cos a sini — V Z: k cos i + Vt,k sin a sin i, (2) 



where i is the inclination of the star (the angle between the 
line of sight and the disc axis) and a is the angle between 
the stellar meridian that contains the line of sight and the 
meridian where the gas parcel under consideration in the k- 
ring is located; each ring has been subdivided into 100 cells 
with Acv = 27r/100 for the calculation of the profile. The 
temperature 71- is taken into account to compute the effect 
of thermal broadening. 

The emissivity of any fc-ring is calculated for the semi- 
forbidden line of Si III] . C III] emissivity was also calculated 
but the density is high and the line is collisionally quenched 
in most of the simulated area thus it was not further stud- 
ied. The atomic parameters for the lines were taken from 
the Chianti Atomic Data Bas^f] (see Appendix A). Colli- 
sional rates have been calculated using thermally averaged 
collision strengths making use of Burgess & Tully (1992) fits 
and rules. For each Tk, the ionisation fraction has been cal- 
culated and the electron density, n e , has been derived. The 
electron temperature, T e , is assumed to be equal to Tk in 
each given fc-ring. 

Si III] emissivity, j, corresponding to the transition be- 
tween levels j and m is given by 



J = 



f(n e ,T e 



"SHII n Si 



(3) 



The population of the relevant levels (rij /nsuu) has been 
computed directly with the Chianti code for each tem- 
perature and electron density. The ionisation fractions 
(nsuu /nsi) have been calculated assuming ionisation equi- 
librium. Solar abundances are used for nsi/nn and the ioni- 
sation fractions to determine uh /n. f(n e ,T e ) is a correction 
factor that takes into account the collisional quenching of 
the linfl 

For a given inclination, the theoretical Si III] profiles 
are computed from the simulation data by adding the con- 
tributions from all the fc-rings with tjj ^ 4R* ; this condition 
was partly set to avoid confusion with any stellar-associated 
component. In addition, the numerical simulations analysed 
(vRB04, vRB06) are not adequate to follow the thermal in- 
stabilities in the boundary layer between the star and the 
disc and, as a consequence, the physical conditions and the 
derived radiative output are not reliable for this region of 
the computational mesh. 

For each simulation in Table 1, the Si III] profile has 
been computed for each time unit (approximately 3 days) 
to track the profile variability, to separate quiescence from 



2 The minus sign has been set to indicate the velocity shift as 
seen by the observer, so negative velocities (or "blueshifts" ) are 
associated with gas motions towards the observer. 

3 URL: www.damtp.cam.ac.uk/user/astro/chianti/ 
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f(n e ,T e ) 



where n e qj tTn and Aj <m are the rates for downward collisional 
and radiative transitions, respectively. Note that when the line is 
collisionally quenched, then /(n c ,T e ) oc n~ 1 T c 1/ ' 2 < T > (T c ). 
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outburst and to compute the line flickering during the qui- 
escence phase. 



4 ANALYSIS OF THE OUTPUT PROFILES 

A given profile is a set of pairs (V ra d, f) where / is the 
total flux radiated by the outflow into the line of sight with 
projected radial velocity in the range [V ra d— 5 km s" 1 , V ra d+ 

5 km s -1 ]. Vrad ranges from -300 km s _1 to 300 km s~\ 
taken in steps of fO km s _1 that oversample the thermal 
broadening (typically ~ 30 km s _1 ) by a factor of 3. 

The system alternates episodic ejections of matter (out- 
bursts) with quiescent periods where accretion is dominant 
(see vRB04). During outburst, the magnetic star-disc link is 
basically broken: the stellar dipolar field lines connecting to 
the disc are opening up and matter is flowing along the open 
stellar and disc field lines forming an outer stellar wind and 
enhanced inner disc wind. In Fig. 2, Si III] emissivity maps 
are shown of the reference model at four different times: the 
first three images are during the quiescence period while the 
last one is during outburst. 

A reliable comparison between the theoretical profiles 
and the observed profiles requires the identification of a set 
of statistically meaningful quantities to describe the profiles. 
We have used for this purpose: 



The velocity of the peak of the line, V p . 



v, J 



• The line asymmetry, A, defined as: A — l b r , where / 

represents the flux at projected radial velocity v, v b marks 
the minimum radial velocity of the blue wing of the line, 
and v r marks the maximum radial velocity of the redwing 
of the line. 

• The line centroid, < V >, defined as: < V >= Jl r , . 

2 2 

• The dispersion, a, defined as: a = %^ f V> ^ . 

Moreover, as the redwing of the profile may have a vari- 
able contribution from mass infall (see RY Tau profiles in 
GdCV07), dispersion and line centroid have also been eval- 
uated only for the blue wing of the profile: 

• The line centroid of the blueshifted emission, < Vb >, 

S° v f*v 

defined as: < Vb >= Jo , . 

• The dispersion in the blueshifted emission, a b , defined 



as: <ib 



In Fig. 3 the temporal variations of these quantities are 
plotted for a whole quiescence-outburst cycle for the "Ref- 
erence" model. The quiescence period can be readily iden- 
tified, as well as the outburst period when the profiles are 
more weighted to the blue and broader than in quiescence. 
In other words, the quantities more sensitive to trace mass 
ejection are: a, o b and < Vb > which, in turn, are rather 
insensitive to the inclination of the source with respect to 
the line of sight. 

The effects in the Si III] profiles due to an outburst are 
shown in Fig. 4 for the reference model: the line broadens 
and the emission region expands towards the boundaries of 
the computational domain. However, during outburst there 
is also some contribution from matter infall (accretion) to 
the profile (see bottom panels in Fig. 4), when some matter 





4v^ 



-200 200 
VJkm/s) 



Figure 4. Theoretical profiles of the Si III] line derived from the 
reference model of the numerical simulations. The time averaged 
profile in quiescence is plotted in the left hand panels and the time 
averaged profile in outburst is plotted in the right hand panels, 
for various inclinations. 



is channelled along the rising and closing field lines into an 
infalling accretion flow - at the same time as enhanced winds 
are starting or ceasing to form. In summary, a clear corre- 
lation between outflow and accretion signatures is predicted 
from these numerical simulations. 

Though these results are extremely promising because 
the correlation between infall and outflow signatures has 
been reported for various cTTSs (i.e. RU Lup by Stempels 
& Piskunov 2002, RW Aur by Gomez de Castro & Verdugo 
2003 and A A Tau by Bouvier et al 2007), they should be 
considered just as a first order approach; the axisymmetry 
assumption imposed in the code and in the calculation of 
the radiative output, forces one to assume that large plasma 
tori are ejected (accreted) along the disc axis at the time of 
outburst (as the smoke rings produced by master smokers) . 
However, in the T Tauri systems mass ejection occurs in 
the form of blobs (plasmoids) from reconnecting magnetic 
fields, which have a much smaller emitting volume than the 
tori assumed in this axisymmetric computation. As a conse- 
quence the fluxes predicted in this paper for outburst (both 
infall and outflow) are overestimated. For this reason, we 
should concentrate on the analysis of the profiles during the 
quiescence phase. 

During quiescence, flickering produces variations of Ob 
and < Vb > that are highly correlated, as shown in Fig. 5 
for two numerical simulations [Reference and Dynamo *) and 
two possible inclinations; the correlation coefficient is be- 
tween 0.96 and 0.9998 depending upon the simulation. This 
high correlation indicates that the shape of the line profile 
does not vary significantly - to the first moment - during 
quiescence, and, in turn, the variations in the profile caused 
by mass ejections are averaged out by the large broaden- 
ing induced by the thermal motions of the gas. The largest 
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Figure 3. Temporal variations of the statistical quantities used to characterise the line profiles in the simulation of the (magnetospheric) 
reference model. The data corresponding to inclinations of 0°, 23°, 45°, 68° and 90° are represented by solid, dotted, short-dashed, 
long-dashed and dot-dashed lines, respectively. The time refers to after the start of the analysed time series. 



variations in at are observed for the pole-on (inclination 0°) 
orientation. 

The time-average values bb and < Vb > are given for 
each simulation and five inclinations during the quiescence 
period in Table 2. The dispersion about these average values 
is also given; note that this dispersion reflects the short-time 
scale variability (the flickering) of the highly non-stationary 
star-disc interaction. The theoretical line flux calculated 
from the simulations is provided in the last column. As no 
extinction effects are taken into account in the calculation of 
the line profiles and fluxes, the fluxes given in the table are 
the same for all the orientations in each modeff). The flicker- 
ing during quiescence produces slight variations in the line 
fluxes in the range of ~ 10%. The fluxes derived are sim- 
ilar to those derived for warm disc models (see Table 2 in 
Paper II). Average profiles for the various models and three 
inclinations are displayed in Fig. 6. 

Several general conclusions can be derived from the 
analysis of the quantities and profiles: 

• Enhancing the dipolar stellar magnetic field enhances 
erf,, < Vt > and their respective dispersions. Thus these 
quantities are good tracers of the magnetic field strength 
and therefore of the wind launching capability. The reason 
is that the centrifugal gear, for a given accretion rate, de- 
pends on the location of the Alfven surface. Stronger fields 
push this surface further away from the star, thereby raising 
the wind terminal speed as well as increasing the radius w 
of the ring radiating in Si III], and thus increasing also the 
dispersion of the projected velocity in the region where the 
emitting gas is located. 

• bb and < Vb > are sensitive to the viewing angle (the 
inclination of the line of sight with respect to the star) 
for strong stellar dipolar fields or for stellar fields that are 
dynamo induced. In general, the larger the inclination the 
larger the values of bb and < Vb >; however, the maxima 
of b b and < Vj, > are achieved at inclination ~ 68°, which 

5 The circumstellar extinction is inclination dependent since it is 
caused by dust lifted from the disc 



Table 2. The 5-j, and < Vb > quantities (averaged over time in 
the quiescence period) derived from the simulations 



Men 


Inclination 


< v b > 


bb 


Flux 

(erg s _1 cm -2 ) 


Reference 


0° 


-64 ± 4 


75 ± 4 


(4.0 ±0.5) x 10~ 8 




23° 


-63 ± 3 


74 ± 3 






45° 


-63 ± 2 


77 ± 3 






68° 


-64 ± 3 


78 ± 3 






90° 


-58 ± 2 


71 ± 2 




Mag-2kG 


0° 


-81 ± 7 


93 ± 8 


(3.9 ±0.6) x 10~ 8 




23° 


-81 ± 6 


96 ± 7 






45° 


-85 ± 7 


104 ± 8 






68° 


-90 ± 8 


107 ± 9 






90° 


-82 ± 8 


97 ± 9 




Mag-5kG 


0° 


-89 ± 14 


103 ± 16 


(3.9 ±0.6) x 10~ 8 




23° 


-90 ± 12 


109 ± 15 






45° 


-101 ± 11 


124 ± 14 






68° 


-108 ± 11 


130 ± 13 






90° 


-101 ± 14 


121 ± 18 




Passive 


0° 


-45 ± 6 


54 ± 6 


(1.9 ±0.2) x 10~ 8 




23° 


-45 ± 6 


57 ± 6 






45° 


-48 ± 6 


63 ± 7 






68° 


-50 ± 6 


65 ± 7 






90° 


-46 ± 4 


60 ± 5 




Dynamo* 


0° 


-64 ± 8 


77 ± 9 


(3.1 ±0.8) x 10~ 8 




23° 


-70 ± 5 


86 ± 6 






45° 


-83 ± 7 


102 ± 8 






68° 


-89 ± 10 


108 ± 11 






90° 


-85 ± 10 


103 ± 11 





roughly corresponds to the angle between the field lines and 
the disc axis in the main ejection region: the current sheet 
region between the disc field and the stellar field. 

• Passive discs produce narrow lines. There are two rea- 
sons for that. As discussed in vRB04, the initial magnetic 
field in the disc due to the penetrating stellar magnetosphere 
is expelled out of the disc and the wind flows nearly parallel 
to the disc surface in the low corona. As a result, the veloc- 
ity dispersion is smaller. In addition, the Alfven surface is 
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Figure 5. Left panels: variation of o~b versus < Vb > and versus 
the line flux (in arbitrary units) during quiescence and inclina- 
tions 68° and 0° (in the top and bottom panels, respectively). The 
temporal line profile variations during quiescence are represented 
in the right panels. Top: "Reference model", bottom: "Dynamo*" 
model. 




-200 200 -200 200 -200 200 



V„,(km/a) V„,(km/s) V,. a (km/s) 

Figure 6. Theoretical profiles of the Si III] line during quiescence 
for the five numerical models analysed (see Table 1) and three 
inclinations with respect to the line of sight: 23°, 45° and 68°. 



5 COMPARISON WITH OBSERVATIONS AND 
SUMMARY 

In Fig. 7, the time-averaged values at and < Vt > are plotted 
for the various models studied in this work. The error bars 
mark the dispersion around the average values in quiescence, 
caused by the flickering. All the inclinations are plotted and 
models are colour coded. The general trends described in 
the previous section are clearly apparent. 

o b and < Vb > are also plotted for warm disc winds 
(purple squares) for the analytical solutions of Gomez de 
Castro & Ferro-Fontan (2005). As expected, warm disc wind 
solutions are less efficient than winds driven from the star- 
disc boundary layer. 

There is a correlation between db and < Vb > , extending 
from the profiles derived in the numerical simulations of this 
work to warm disc winds (to some extent). This correlation 
traces the increasing efficiency of MHD centrifugally driven 
winds and shows that larger terminal speeds are associated 
with larger bulk motions at the base. Notice that inclination 
effects are difficult to sort out in this diagram. 



located at the innermost boundary of the star-disc system 
and thus the centrifugal gear is negligible. 

• Active stellar dynamos are able to generate a rapid out- 
flow even for a moderately strong magnetic field. The rea- 
son is because in the "Dynamo *" model the stellar wind is 
mostly driven by both poloidal and toroidal magnetic pres- 
sure, as well as by gas pressure. 



5.1 The observed Si III] profiles 

The information available on the Si III] profiles of TTSs 
come from the observations carried out by Gomez de Castro 
ct al (2003, hereafter GdCVFF03) with the HST. The Si III] 
profiles are plotted in Fig. 8, as well as the C III] lines needed 
for the density diagnostic. All the targets (unless AK Sco) 
have reported jets allowing a clear identification of the out- 
flow signature (see Table 3) . As the size of the spectrograph 
slit was 0.2 arcseconds, any extended structure within 28 
AUs could contribute to the line flux. 

The contribution from the large scale jet to the line flux 
can be discriminated by the density of the emitting plasma. 
The electron density in the jet decreases very rapidly with z: 
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Table 3. Properties" of the TTSs observed in Si III] 




Figure 7. Representation of cr;, versus < > for all models 
and inclinations of Table 2. The time-averaged values erf, and 
< Vj > are represented as squares, while the error bars indi- 
cate true temporal variations of the profiles during quiescence. 
Dark blue, black, green and light blue are used to represent the 
models: Passive, Reference, Mag-2kG and Mag-5kG, respectively. 
Red is used to represent the Dynamo* model. Purple squares at 
the bottom right represent predictions from the warm disc wind 
models by Gomez de Castro & Ferro-Fontan (2005). The observed 
values for the TTS: RY Tau are plotted as two big black triangles, 
representing the two extreme values that can be derived from the 
observations carried out with the HST in 1998 (Gomez de Castro 
& Verdugo 2001, 2007). 



from ~ 10 9 - 10 10 cm -3 at the base, to ~ 10 4 - 10 6 cm" 3 in 
the jet due to the rapid radial expansion of MHD centrifu- 
gally driven winds. At the low densities of the jet, the forbid- 
den lines of [Si III](A1983) and [C III](A1907) are stronger 
than their semiforbidden counterparts. Also, the ratio Si- 
IIIJ/CIII] is < 1 (Keenan et al 1992), as observed in the 
UV spectra of protostellar jets and Herbig-Haro objects (see 
Gomez de Castro & Robles 1999 for a compilation). 

This is not observed in the HST/STIS spectra of TTSs 
as shown in Fig. 8 (from GdCVFF03). Typically, the Si III] 
and C III] lines have similar strengths. The forbidden [C III] 
line is only detected in T Tau and RU Lup and, only in 
T Tau, has a non negligible contribution to the Si III] flux. 
Thus Si III] emission is produced at the base of the outflow 
in most of the sources. 

The Si III] flux may include additional contribu- 
tions from accretion shocks and the stellar magneto- 
sphere/atmosphere. Discriminating between these three 
components is much more difficult since they share similar 
temperature and density regimes. 

RY Tau is the only TTS that has been observed twice 
in the Si III] lines allowing to separate cleanly the contri- 
bution from accretion and outflow becoming an ideal test 
for the models. a b and < V b > for RY Tau (from GdCV07) 
are also plotted in Fig. 7. Notice that the line centroid is 
too blueshifted for the observed dispersion. This is because 
the contribution from the outflow is dominated by a narrow 
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References: (1) Hartigan et al 1995; (2) Hirth et al 1997; 
(3) Takami et al 2001; (4) Clark & Bouvier 2000; 
(5) Andersen et al 1989; (6) Ardila et al 2002; 
(7) Stempels et al 2007; (8) Muzzerolle e al 2003. 
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Figure 8. Measured Si III], C III] and [C III] profiles of TTSs 
(from GdCVFF03). The rest velocity corresponds to the rest 
wavelength of the Si III] and C III] lines in the left and right 
panels, respectively. The "velocities" of the [C III] and C III] 
lines are marked with dashed lines in the right panel for clarity. 



blueshifted component with a long tail to redwards shifted 
velocities (see Fig 1 in GdCV07) . This disagreement between 
theory and observations can be caused by the lack of cylin- 
drical symmetry of the outflow (i.e. mass ejected as blobs 
from the current layer). 



5.2 The a versus < V > 

The a versus < V > diagram is better suited to break the in- 
clination degeneracy and examine the relevance of outflows 
with respect to magnetospheric processes in TTSs. The cen- 
troid of the line emission from the wind moves to higher 
blueshifts as the inclination with respect to the line of sight 
decreases however, no significant inclination effects are ex- 
pected from magnetospheric models. 

The values of a and < V > derived from the Si III] 
profiles (see Fig. 8) are plotted for comparison. Notice that 
< V >~ for most sources and that there is not any trend 
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related with the inclination. This suggests that that Si III] 
emission is dominated by circumstellar structures rather 
than by disc winds. In fact, RW Aur, the star with the 
highest dispersion, is surrounded by a torus/belt of ionized 
plasma (Gomez de Castro & Verdugo, 2003). However, even 
stars that have not this kind of circumstellar structures, as 
DE Tau, display profiles with large broadenings that can- 
not be accounted solely by stellar rotation. Thus, it seems 
that the dominant contribution to the Si III] emission is 
produced in extended magnetospheric structures where the 
large broadening traces either the accretion flow (Hartmann 
2009) or with large macro-turbulence fields, as the hypoth- 
esised to reproduce the broad profiles of AK Sco (Gomez de 
Castro 2009). 

A small contribution from the outflow to the Si III] flux 
is detected RU Lup, T Tau and RY Tau (once the contribu- 
tion from accretion is subtracted). It shows as a bluewards 
shift of the line centroid. Unfortunately, further observations 
are required to disentangle magnetospheric and wind contri- 
butions prior to an accurate diagnostic of the wind drives. 

We would like to conclude this work with a word of cau- 
tion. This research is limited by the characteristics of the 
simulations analysed, such as the axisymmetry, the rather 
small extent of the computational domain and the very 
efficient tapping of thermal energy imposed by the poly- 
tropic approach in the simulations. More realistic models 
implementing radiative cooling and adaptive mesh refine- 
ment techniques to track the thermal instabilities are cur- 
rently under development. These new simulations will pro- 
duce more detailed models of the dynamics and the radia- 
tive output from the disc-star interface; however, they will 
not modify the strong dependence of the profiles on the in- 
clination shown in Fig. 8, which is caused by the coupling 
between ejection and rotation and holds even on scales of 
0.1 AU. 



i=90° 

i = 6B" 




-100 -80 -60 -40 -20 20 
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Figure 9. Representation of a versus < V > for all models and 
inclinations of Table 2. The time-averaged values a and < V > are 
represented as squares, while the error bars indicate true tempo- 
ral variations of the profiles during quiescence. Dark blue, black, 
green and light blue are used to represent the models: Passive, 
Reference, Mag-2kG and Mag-5kG, respectively. Red is used to 
represent the Dynamo* model. Purple squares at the bottom right 
represent predictions from the warm disc wind models by Gomez 
de Castro & Ferro-Fontan (2005). The observed values for all the 
TTSs observed with the HST/STIS (DE Tau, AK Sco, RY Tau, 
RW Aur, T Tau and RU Lup) are plotted as big black triangles. 
The two observations of RY Tau are plotted; the less blueshifted 
include the contribution from the accretion flow as in the observ- 
ing campaign in 2001. 
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APPENDIX A: ATOMIC PARAMETERS FOR 
THE CALCULATION OF THE LINE 
EMISSIVITIES 
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Table Al. Atomic parameters of the semiforbiddcn lines 



Element 


Wavelength 


3 


m 


A ■ 




Type« 


C< 2 > 








(3) 












s- 1 








0.0 


0.25 


0.50 


0.75 


1.00 


Si III 


1892.033 


3 


1 


1.544X10 4 


3 


2 


0.4 


2.931 


1.655 


1.063 


0.4002 





C III 


1908.737 


3 


1 


97.32 


3 


2 


1.1 


0.3688 


0.3419 


0.3047 


0.2102 


0.0888 



Transition type according to Burgess & Tully (1992) classification. Type 2 includes transitions which are induced by 
either an electric multipolc or a magnetic multipole, e.g. forbidden transitions. 

Collisional rates have been computed following Burgess & Tully (1992). Thermally averaged collision strengths, 
< ~"~3,m > (x), are given by interpolation for each temperature, T c , using the coefficients indicated in the table with 
x(T c ) = (k B T e /E jtm )/(C + (k B T e /E jtm )) and E jiTn = hc/X j:Tn . The collision rate is: 8.6233 X 10" 6 n e T e _1/2 < T jiTt 
(ks- Boltzman constant; h: Planck constant; c: speed of light; wj: statistical weight of level j.) 



